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Abstract 



We have measured the polarization of J /if) and ip{2S) mesons produced in 

pp coUisions at ^/s = 1.8 TeV, using data collected at CDF during 1992-95. 

The polarization of promptly produced J/^p [^'(25')] mesons is isolated from 

those produced in B-hadron decay, and measured over the kinematic range 

4[5.5] < Pt < 20 GeV/c and \y\ < 0.6. For Pt > 12 GeV/c we do not observe 

significant polarization in the prompt component. 
PACS numbers: 13.85.Qk, 13.88.+e, 13.20.Gd 
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The production of heavy quarkonia states, cc and bb, provides a useful system for the 
study of Quantum Chromo dynamics (QCD), as it involves both perturbative and nonpertur- 
bative energy scales. In pp collisions, charmonium production occurs through three mech- 
anisms: direct production, the decay of heavier charmonia, and the decay of 6-fiavored 
hadrons. The first two mechanisms are collectively known as "prompt" because they are 
observed to occur at the pp interaction point. 

The Collider Detector at Fermilab (CDF) collaboration previously reported results on 
the production of J/ip and ipi'^S) mesons [0,11 . The measured cross sections for direct 
production were on the order of 50 times larger than predicted by the Color Singlet Model 
(CSM) 0. However, calculations based on the Nonrelativistic QCD (NRQCD) factorization 
formalism [^,^ are able to account for the observed cross sections by including color octet 
production mechanisms. This leads to the prediction that directly produced ip mesons will 
be increasingly transversely polarized at high Pt (In this Letter we use to denote 

either J/tp or ip{2S) mesons.) This is because the production of ip mesons with Pj- ^ is 
dominated by gluon fragmentation. It is predicted that the gluon's transverse polarization 
is preserved as the cc pair evolves into a bound state ip meson. On the other hand, the 
Color Evaporation Model predicts an absence of polarization In this Letter, we report 
on measurements of the polarization of promptly produced ip mesons at CDF. Our analysis 
also yields as a byproduct the effective polarization of the ip mesons produced in i?-hadron 
decays. 

CDF is a multi-purpose detector designed to study high energy pp collisions produced 
by the Tevatron ||^. The CDF coordinate system is defined with the z axis along the 
proton beam direction. The polar angle 6 is defined relative to the z axis, r is the per- 
pendicular radius from this axis, and (p is the azimuthal angle. Pseudorapidity is defined as 
rj = — ln[tan(^^/2)]. Three charged-particle tracking detectors immersed in a 1.4 T solenoidal 
magnetic field surround the beamline. This tracking system is contained within a calorime- 
ter, while drift chambers outside the calorimeter identify muon candidates. 

The innermost tracking device is a four-layer silicon microstrip detector (SVX) located 
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at radii between 2.9 and 7.9 cm from the beam axis. The SVX is surrounded by a set of 
time projection chambers (VTX) extending out to a radius of 22 cm. An 84 layer cyhndrical 
drift chamber (CTC) measures the particle trajectories in the region between 30 and 132 cm 
from the beam. This tracking system has high efficiency for detecting charged particles with 
momentum transverse to the beam Pt > 0.4 GeV/c and \r]\ < 1.1. Together, the CTC and 
SVX measure charged particle transverse momenta with a precision of (Tp^/Pt = 0.007© 
0.001 • Pt (with Pt in GeV/c). The impact parameter resolution is Ud — (13 + 40/Pt) A*ni 
for tracks with SVX and CTC information. 

The central muon detection system consists of four layers of planar drift chambers sepa- 
rated from the interaction point by five interaction lengths of material. This system detects 
muons with Pt ^ 1.4 GeV/c and |?7| ^ 0.6. Dimuon candidates used in this analysis are 
collected using a three-level trigger. The first level trigger requires that two candidates 

be observed in the muon chambers. For each muon candidate the first level trigger efficiency 
rises from ~40% at Pt = 1-5 GeV/c to ~93% for muons with > 3.0 GeV/c. The second 
level trigger requires one or more charged particle tracks in the CTC, reconstructed using 
the central fast track processor (CFT). The CFT performs a partial reconstruction of all 
charged tracks with Pt above ~2 GeV/ c. Muon candidates found by the first level trigger 
are required to match a CFT track within 15 degrees in azimuth. The third level trigger per- 
forms three-dimensional track reconstruction and accepts dimuon masses in a broad window 
around the J/'4> and 11^(23) masses. 

The data used in this study correspond to an integrated luminosity of 110 pb^^ and 
were collected between 1992 and 1995. Following the online data collection, additional 
requirements are made offline to identify the signals and to reduce the backgrounds. To 
identify muon candidates and reduce the rate from sources such as n/K meson decay-in- 
flight, we require that each track observed in the muon chambers be associated with a 
matching CTC track. These matches are required to pass a maximum cut of 9 and 12 
(for 1 degree of freedom) in the and z views respectively. Also, we require Pt greater than 
about 2 GeV/c for each muon candidate. This requirement ensures that the muon trigger 
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and reconstruction efficiencies are well understood, to avoid biases in the decay angular 
distributions of the charmonia states studied below. 

The measurement of the polarization of ip mesons is made by analyzing their decays to 
fi~^fi~ in the helicity basis, in which the spin quantization axis lies along the ip direction in 
the pp center-of-mass (lab) frame. We define 6* as the angle between the /i"*" direction in 
the ip rest frame and the ip direction in the lab frame. The normalized angular distribution 
I {cos 6*) is given by 

I{cose*) = ^^^^{l + a cos'e*) (1) 

Unpolarized ip mesons have a = whereas a = +1 or —1 correspond to fully transverse or 
longitudinal polarizations respectively. Experimentally, the acceptance is severely reduced 
as I cos 6** I approaches 1, due to the Pt cuts on the muons. Our method for determining a is 
to fit the observed distributions of cos 6* to distributions derived from simulated ip fi'^fi~ 
decays. The Monte Carlo simulation accounts for the geometric and kinematic acceptance 
of the detector as well as the reconstruction efficiency as a function of cos^*. 

In order to extract the polarization parameter a for promptly produced ip mesons, we 
separate the prompt component from the 5-decay component using the proper decay length 
of each event. For ip candidates with one or both muons reconstructed in the SVX (the 
SVX sample), we define a vector pointing from the pp collision point to the tp decay vertex. 
The transverse decay length Lxy is then defined as the projection of this vector onto the 
ip transverse momentum. The proper decay length ct is related to the transverse decay 
length by ct = {M^L,y)/{F^^„P^), where is the i; mass. Here F^^^^ is a correction 
factor obtained from Monte Carlo studies [jlO|, which accounts for the fact that we are 



using the ip Pt instead of the i?-hadron P^. Prompt events have ct consistent with zero 
whereas 5-decays have an exponential ct distribution; the detector resolution smears the ct 
distribution. We fit the ct distribution to obtain the relative fractions of prompt and -B-decay 
production. Details of this ffiting procedure are given in The measured fraction of J/ip 
mesons which come from S-hadron decay increases from (13.0 ± 0.3)% at P^^"^ = 4 GeV/c 
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to (40 ± 2)% at 20 GeV/c. For ip{2S) mesons, an increase from (21 ± 2)% to (35 ± 4)% is 
seen in the range from 5.5 to 20 GeV/c. 

The proper decay length measurement allows us to divide the data into two samples: a 
short-lived sample dominated by prompt production, and a long-lived sample dominated by 
-B-decays. The short-lived sample is defined by —0.1 < ct < 0.013[0.01] cm, and the long- 
lived sample by 0.013[0.01] < ct < 0.3 cm, for the J/ip [^'(25')] analysis respectively. The 
boundary between the two ct regions has been optimized separately for the J/ip and ip{2S) 
samples, to maximize the purity of prompt decays in the short-lived sample and -B-decays 
in the long-lived sample. Depending on Pj^, the prompt fraction in the short-lived sample 
ranges from 85% [86%] to 96% [95%], and the 5-decay fraction in the long-lived sample ranges 
from 83%[86%] to 98%[91%], for J/ip['ip{2S)] respectively. 

The J/ip polarization is measured in seven bins, covering a range of 4 — 20 GeV/c. 
Using a 3 standard deviation mass window around the J/ ip peak, our data sample consists of 
180,000 signal J/ip events, with a signal to background ratio of about 13. The J/ip sample is 
divided into three subsamples: the short-lived and long-lived SVX samples described above, 
and a third sample (the CTC sample) in which neither muon has SVX information and no 
ct measurement is made. In each P^^^ bin, we measure the prompt polarization (ap) and 
the effective polarization of J/ip mesons from i?-hadron decays {as)- (We refer to as 
"effective" because 6* is defined using the lab frame, not the i?-hadron rest frame — in effect 
this dilutes any polarization from the S-decay toward zero.) We find that it is not feasible 
to make separate polarization measurements for direct J/ip production and for production 
from Xc and ip{2S) decays. The latter sources account for 36±6% of the prompt component, 
with only a small P^^^ dependence [0. 

The J / ip polarization is measured by fitting cos 9* distributions in data to a set of Monte 
Carlo templates |[11[]. The templates are generated by processing simulated samples of 
J/ip ^ fj^~^fi~ decays with a detector and trigger simulation. The polarization is obtained 
using a fit of the data to a weighted sum of transversely polarized and longitudinally 
polarized templates. The fitted weights yield the polarization. Two transverse/longitudinal 
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template pairs are generated, using measured prompt and 5-decay spectra f^. The 
cos 6** distribution of background events is modeled in the fit using sidebands around the 
J/xf} mass peak. The fit is performed simultaneously on the SVX short-lived, SVX long- 
lived and CTC samples, with two fit parameters: ap and as- To account for the mixture 
of prompt and -B-decay components in each sample, the relative fractions of prompt and B- 
decay templates in each are fixed in the fit using the results of the lifetime fit. The B-decay 
fraction in the CTC sample is assumed to be the same as in the SVX sample, because the 
two samples differ primarily in the z position of the primary vertex. Within each P^^^ bin, 
a small correction is applied to the P^^^ distributions of the Monte Carlo samples so that 
they match with those in the data. As an example, the fit in the Pt range 12-15 GeV/c is 
shown in Fig. |l|. 

Three sources of systematic uncertainty are evaluated: the trigger efficiency, the fitted 
prompt and -B-decay fractions, and the P^^^ spectra used in making the Monte Carlo tem- 
plates. Except in the lowest Pt bins, the systematic uncertainties are much smaller than 
the statistical uncertainties. Our fit results are listed in Table |, and ap is compared with 
a theoretical NRQCD prediction in Fig. |^. 

The measurement of the ipi'^S) polarization is made in three Pt bins covering 5.5 — 
20.0 GeV/c. Both muons are required to be reconstructed in the SVX. The resulting dimuon 
mass distribution is fitted with a Gaussian signal and a linear background. We find a total 
of 1855 ±65 signal 4'{2S) events, with a signal to background ratio of about 1 in a 3 standard 
deviation mass window around the ip{2S) mass. 

As discussed above, the sample in each Pt bin is further divided into two subsamples 
based on the ct distribution. Because the statistics are lower than in the J/ip case, we use 
10 bins in | cos 6* \ . The number of signal events in each | cos 6* \ bin is obtained by fitting 
its mass distribution. The resulting | cos 9* \ distributions in the two ct subsamples are fitted 
simultaneously to the predicted number of events to extract the ip{2S) polarizations for 
prompt and -B-decay production. The number of predicted events in each | cos 9* \ bin is 
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derived by weighting the normahzed angular distribution J(cos^*) with the detector accep- 
tance [|T2|. We use the measured prompt and i?-decay P^^'^^^ distributions to calculate 
the acceptance. As in the J/ip case, there is a small correlation between the measured 
p^(2S') jjg^j^j]3^^jQ^g QT^(^ ^i^Q polarization. A correction is applied iteratively in the fits to 
account for this dependence. Figure ^ shows the observed angular distributions with their 
polarization fits for the short-lived sample in the three P^^"^^^ bins. 

Three sources of systematic uncertainty are considered: the uncertainty in the event yield 
from the mass fits in the | cos 6* \ bins, the uncertainty due to the error on the fitted prompt 
and i?-decay fractions, and the uncertainty on the | cos^*| acceptance from the Monte Carlo 
modeling of the P-^^^'^^ distributions. The uncertainty due to the trigger efficiency is neg- 
ligible in the Pf^^^^ range used. The systematic uncertainties are much smaller than the 
statistical uncertainties. The fitted values of ap and as a function of P^^"^^^ are listed in 
Table and ap is shown in Fig. || with the NRQCD predictions f^J^l- 

In conclusion, we have measured the polarization of J/tp and ip{2S) mesons produced in 
1.8 TeV pp collisions. The polarization from i?-decays is generally consistent with zero, as 
expected. In both the J/ip and ipi'^S) cases, we do not observe increasing prompt transverse 
polarization at Pt ^ 12 GeV/c. Our measurements are limited by statistics, especially 
for the ip{2S), but they appear to indicate that no large transverse prompt polarization is 
present at high Pj-, in disagreement with NRQCD factorization predictions. 
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TABLES 



TABLE I. Fit results for J/ij; polarization, with statistical and systematic uncertainties. 



Pt bin (GeV/c) 


Mean Pt (GeV/c) 


ap 


ap 


4- 5 

5- 6 

6- 8 
8-10 
10-12 
12-15 
15 - 20 


4.5 

5.5 

6.9 

8.8 

10.8 

13.2 

16.7 


0.30 ±0.12 ±0.12 
0.01 ±0.10 ±0.07 

0.178 ±0.072 ±0.036 
0.323 ± 0.094 ± 0.019 
0.26 ±0.14 ±0.02 
0.11 ±0.17 ±0.01 

-0.29 ± 0.23 ± 0.03 


-0.49 ±0.41 ±0.13 
-0.18 ±0.33 ±0.07 
0.10 ±0.20 ±0.04 
-0.06 ±0.20 ±0.02 
-0.19 ±0.23 ±0.02 

0.11 ±[]:i ±0.02 

-0.16±[]:i| ±0.05 


TABLE IL Fit results for ip{2S) polarization, with statistical and systematic uncertainties. 


Pt bin (GeV/c) 


Mean Pt (GeV/c) 


ap 


ap 


5.5-7.0 
7.0 - 9.0 

9.0 - 20.0 


6.2 
7.9 

11.6 


-0.08 ± 0.63 ± 0.02 
0.50 ±0.76 ±0.04 

-0.54 ±0.18 ±0.04 


-0.26 ± 1.26 ±0.04 
-1.68 ±0.55 ±0.12 

0.27 ± 0.81 ± 0.06 
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FIGURES 
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FIG. 1. The J/ijj polarization fit to cos^* distributions in the 12-15 GeV/c bin. Points: side- 
band-subtracted data in the SVX short-hved, SVX long-lived and CTC samples. Dashed lines: 
fit. 
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FIG. 2. (a) The fitted polarization of prompt J/ij: mesons for ly^/'^l < 0.6. Full error bars 
denote statistical and systematic uncertainties added in quadrature; ticks denote statistical er- 
rors alone. The shaded band shows an NRQCD factorization prediction which includes the 
contribution from Xc and ip{2S) decays, (b) The fitted polarization of prompt ip{2S) mesons for 



.,V'(25)| 



< 0.6. Error bars denote statistical and systematic uncertainties added in quadrature. 



Shaded bands show two NRQCD factorization predictions [^,|7|. 
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FIG. 3. Fits to I cos 6*\ distributions in the short-lived i>{2S) data sample, in the three Pr bins. 
Points: data, dashed lines: fit. The acceptance extends further out in | cos6*\ as Pr increases. 
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